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Ghost imaging capability was initially attributed only to entangled photon pairs, ^ but in 
more recent works it has been ascertained also for classically correlated fields, both in the 
single photon regime^' ^ and in the continuous variable regime.^ The main result of all these 
schemes is the non-local imaging of objects both in near- and far-field.^ A general ghost- 
imaging scheme involves a source of correlated bipartite field and two propagation arms 
usually called Test (T) and Reference (R). In the T-arm the object to be imaged is inserted 
and a bucket (or a pointlike) detector measures the light transmitted by the object. The 
R-arm contains an optical setup suitable for reconstructing the image of the object (or its 
Fourier transform) and a position-sensitive detector. 
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In this Letter we show a ghost imaging experiment performed by using the intense fields 
generated by parametric downconversion (PDC) seeded with multimode chaotic hght. As 
usual in ghost imaging experiments, the technique adopted for the retrieval of the image 
consists in the evaluation of the fourth-order correlation between the fields at the detection 
planes. 

In Fig. [T]we show our experimental setup. We injected a coherent coUimated pump beam 
{Xp = 532 nm) and a seed beam (A = 1064 nm) into a crystal of /5-BaB204 (BBO, cut angle 
22.8 deg, 10 mmxlO mmx3 mm, Fujian Castech Crystals) in type-I interaction geometry. 
Both pump and seed beams were provided by an amplified Q-switched Nd:YAG laser (GCR- 
4, 10 Hz repetition rate, 7 ns pulse duration of the fundamental pulse, Spectra-Physics). The 
seed field was randomized by passing it through two independently rotating ground glass 
plates,^ Pi and P2 in order to obtain a pseudo-thermal statistics. The speckles on the object 
plane resulted to be about 320 /im in diameter as evaluated by spatial autocorrelation. The 
object was a hole of 1.6 mm diameter crossed by a straight wire of 0.5 mm caliber. As the Dt 
and D/j detectors, we used a single CCD camera (CA-D1-256T, 16 yumxl6 /xm pixel area, 
12 bit resolution, Dalsa). The CCD sensor, also depicted in Fig. [U for one half recorded 
the single-shot intensity Maps in the R-arm, X/j(x/j), and for the other half realized the 
bucket detector of the T-arm. In the following, X^ designates the sum of the contents of 
110 X 110 pixels of Dy in a single shot. The T-arm includes free propagation over di = 60 cm 
from BBO to object, transmission function t(x), and light collection optics in front of Dt- 
The R-arm contains free propagation over ^2 = 60 cm from BBO to lens L (/ = 400 mm) 
and from L to the CCD camera [q = 60 cm). Distances di, ^2 and q and focal length / 
must satisfy l/{di + ^2) + l/q = 1// to give an imaging system with magnification factor 
M = q/{di + d2)} 

The evolution of the system is described by the unitary operator S generated by the 
multimode PDC hamiltonian^ and the input-output relations that link the output operators 
6i_q to the input operators Oj^q are 

fc.,q = S^a.,^S = U,,^a.,^ + K,qe'^<>4__q z = R,T; t y^ j, (1) 

where Uf^^ — V^^ = 1 and {U, V, y^jq = {U, V, (p}^^- We take a/j q in the vacuum state and 
aT,q in the multimode thermal state represented by the density matrix 

Pr = n ^ IZ Pr,q(m) \m) (m|^q \ , (2) 

q 1 771=0 ) 

where I'm) is the Fock state with m photons in mode q and pT,q{fn) = n^ /{I + nth,q)"^~^^ 
is the photon number distribution on the single mode for nj/i q mean photon number. The 
the output state is obtained by exploiting the Baker- Campbell-Hausdorff formula^' ^ 

Pout = S{pT^\0){0\j,)S^ 




Fig. 1. Experimental setup: Pi, 2? independently rotating ground glass plates: 
Pol, polarizer; t(x), object transmission function; L, lens; tuning angle 24.7 
deg. 
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'T'PDCq = ^,q being the mean photon number per mode generated by spontaneous PDC. 
The marginal distributions on R- and T-arms are multithermal. 

By direct application of the Peres- Horodecki- Simon criterion^"^^ it can be demonstrated that 
the state pout is inseparable for any value of nth,ci-^^ 



The reconstruction of the image is achieved by the computation of the correlation function, 
G'''^^(x/j). This procedure is equivalent first to evaluate the correlation function between 
Ir{xr) and It{^t) 



g^'\^R,^T) = {Ir{^r)It{^t)) - {Ir{^r)){It{^t)), 



(5) 



where (/j(xj)) = (q(xj)cj(xj)) {i = R,T) is the mean intensity of the i-th. beam and 
(/r(xr)/t(xt)) = (4(xfl)cR(xR)4(xr)cr(xT)) (where (...) = Tr[...pr® |0) (0|^]), and then 
to perform the integration 



G'^'\^R) = Jd^Tg'-'\^R,^T). 



(6) 



The propagation to Dr and Dt is described by the corresponding impulse response functions 
hR{'XR,x.'j^) and /it'(xt,x^), whose derivation is straightforward/'^ and the field operators at 
the detection planes become Cj(xj) = /dx^/ij(xj,x96j(x^), where 6j(x) oc Z]qexp(2q ■ x)6jq. 
The factorization rule for (6|j(x'^)6r(x'^)67.(x^)6t(x^)) in this case is exactly the same as 
that for spontaneous PDC/^ thus in the degenerate case Eq. ([5]) becomes 
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If the coherence area of the multithermal field is much smaller than the object, we can 
proceed as in Ref.^^ and get 
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In our experiment we estimate ripoc — 0.5 and nth — 10^^ per spatial mode. The correlation 
function in Eq. (^ reproduces the object and depends on the properties of our source through 
the mean photon numbers of downconverted fields and of the chaotic seeding field. 

The experimental results are shown in Fig. [2j panel a) displays a typical sample of chaotic 
intensity Map XR(xi^) recorded in single shot, whereas panel b) contains the map of the 
correlation coefficients 
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evaluated over an ensemble of 9500 single shot Maps, where a^(/) = (P) — (J)^ and both 
intensities and variances were evaluated by subtracting the statistical contributions of the 
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Fig. 2. a) Single-shot Map recorded in the R-arm; b) ghost image; c) contour 
plot of the local visibility. 



background noise. To measure the noise we considered the values recorded by the CCD in 
a non-illuminated sensor region. Panel b) in Fig. [2] shows that the recovered image has the 
correct size, as compared to the original object, being M = 0.5. 

To quantify the quality of our technique, we define the local visibility of the reconstructed 
image (see Ref.^^) as 
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{Ir{^r)It) (/ii(x^))(/r)+G(2)(x^) 

and evaluate it for the ghost image in Fig. [2]b). Figure [2] c) shows the contour plot of V(x/j) 
at the marked levels including the maximum value 0.046. 

A tradeoff between visibility and resolution of the reconstructed image has been reported^^ 
that depends on the number of spatial modes (coherence areas, Ac) that illuminate an object 
of area Aq- In our case Aq/Ac — 1600^/320^ = 25. This number being rather low, we observe 
a quite poor image resolution and a quite high visibility. To confirm our results we calculate 
the number of spatial and temporal modes, involved in the interaction. We start by studying 
the statistical distributions of the beam intensities detected either in a single pixel x^j or by 
the bucket detector in the T-arm. In Fig. [3] we show the experimental distributions along with 
the best fitting curves obtained by convolving the theoretical multithermal distribution^' ^^ 
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with the experimental one obtained for the background. Here X indicates the single-shot 
values of either Xr(xr) or It, whose mean values and variances obviously coincide with 
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Fig. 3. a) Statistical distribution of the intensity in a single pixel of the D^ 
detector (dots) and multithermal fit giving fiji = 1 (full line); b) statistical 
distribution of the intensity on the bucket detector (dots) and multithermal 
fit giving fiT = 35 (full line). 



those of the operators J/?(x/j) and It- From the data in panel a) we obtained /ir = 1, 
independently of the position x/j we tested, and from the data in panel b), /i^ = 35. As 
the statistical distribution in Fig. [3] a) reflects the behavior of the field from shot to shot, 
we interpret fiR as the number of temporal modes in the field. ^'^^ On the other hand, we 
can interpret fiT as the product of the number of temporal modes times that of the spatial 
modes in the area covered by the 110 x 110 pixels of D^- As /i/j ~ 1 we conclude that we 
have a number of spatial modes in the bucket, 35, greater than the number, 25, of those 
illuminating the object (see above): this is to be expected as the bucket integration covered 
an area wider than Aq- From Eq. ( ITTl) we get (t'^{Ir^t) = {lR,T)'^/fJ'R,T, which allows linking 
the visibility of Eq. fITOl) to the number of thermal modes in the incoherent T- and R-beams. 
By using Eq. ([9]) we can rewrite Eq. flTOj) as 
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and make an independent estimation of ^/JTbPt as a function of the experimental values of 
correlation coefficients and visibility. By this method we find ^/JTrJIt — 6.8. As from the fits 
in Fig. Owe found the value ~ 5.9, the experimental results are self-consistent. 
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In conclusion, we have implemented a new source of correlated beams suitable to perform 
ghost imaging and, as we expect, ghost diffraction experiments. Besides the advantage of 
using detectors that measure intense light, the major benefit of this source as compared to 
those operating on spontaneous PDC is the possibility of tuning resolution and visibility 
of the ghost image, for instance by modifying the spatial coherence properties of the seed 
beam. In fact we observed that the spatial-coherence structure of the seed is preserved on 
both arms of the seeded PDC emission because all its spatial components undergo interaction 
with the pump. The only limitation could arise from the competition between the speckle 
divergence (< 0.3 mrad) and the angular bandwidth of the interaction (~ 5 mrad). Note that 
the contribution of spontaneous PDC, which is expected to have a greater angular spread, 
is negligible in our experiment. 

We thank M.G.A. Paris (University of Milan) for fruitful discussions. 
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